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Abstract 

The fundamental tenet of the classical supernovae-driven wind model of elliptical 
galaxies is that the residual thermal energy of all supernovae remnants (SNRs) provide 
sufficient energy to overcome the binding energy of the remaining interstellar gas, thereby 
driving a global galactic wind. We re-examine model predictions of this epoch of wind ejec- 
tion ^GW, highlighting a heretofore underappreciated sensitivity to the adopted remnant 
thermal energy formalism, and illustrating cases in which previous work may have sub- 
stantially overestimated ^gw- Arguments based upon chemical evolution alone, put forth 
to reject the hypothesis of dark matter distributions similar to the luminous component 
in spheroids, are shown to be tenuous. Finally, the predicted enrichment of intracluster 
gas during the wind phase of cluster ellipticals, and its relation to the selected SNR in- 
terior thermal energy evolutionary scheme, is addressed. Despite the success of previous 
wind models, our results still call into question the correctness of the simple analytical 
approach used thus far, and imply that a more appropriate technique should be adopted 
in the future. 
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1. Introduction. It has long been recognised that the colours of elliptical galaxies 
become increasingly red with decreasing magnitudes, or conversely, increasing mass. This 
effect has been commonly attributed to an increasing average galactic metal content with 
increasing galactic mass (e.g. Faber 1977, for an early review). The key to understanding 
the origin of these correlations was provided by Mathews & Baker (I97I), but not fully 
appreciated until Larson (1974). 

In order to explain the paucity of interstellar gas in most ellipticals Mathews & 
Baker (1971) postulated that much of it had been strongly heated by supernovae (SNe) 
explosions and driven out by a hot galactic wind once the residual thermal energy of all 
SNe remnants exceeded the binding energy of the remaining gas, thereby bringing to a halt 
the bulk of active star formation. Subsequent evolution would then be regulated by the gas 
returned to the interstellar medium (ISM) from dying stars. It was left to Larson (1974) 
to recognise that the binding energy per unit mass of gas is higher in the more massive 
galaxies, and thus these systems would retain their gas for longer initial periods before 
reaching this epoch of galactic wind ^gw, thereby attaining higher metallicities in a manner 
consistent with the observed mass-metallicity relationship. Further developments over the 
past decade can be attributed to Saito (1979), Arimoto & Yoshii (1987), Matteucci & 
Tornambe (1987), Angeletti k Giannone (1990), Bressan, Chiosi k Fagotto (1994), Elbaz, 
Arnaud k Vangioni-Flam (1994), and Gibson (1995). It is generally accepted that wind 
models of this nature may be inadequate in some details, but are still the framework of 
choice for modeling the early evolution of spheroids (Renzini et al. 1993). 

Determining the epoch of global gas ejection is an imperative aspect of these wind 
models as ^gw governs the cessation of the bulk of star formation, as well as the amount 
and metallicity of expelled gas, and the current spectrophotometric properties of ellipticals. 
Angeletti k Giannone (1991) have discussed the sensitivity of ^gw predictions to the 
assumed initial proto-galactic radius (contributing to the binding energy of the spheroid), 
and Matteucci (1994) has explored recently the importance of the astration parameter v 
(contributing to the timescale for star formation). Our plan herein is to re-examine briefly 
one of the fundamental ingredients to all supernovae-powered wind models, speciflcally, 
the formalism adopted to represent the evolution of a supernova remnant's (SNR) internal 
thermal energy Stiiit), and the availability of such energy for driving a galactic wind. 

In Section 2 we outline the basic chemical evolution and wind model framework used 
in our study. We present the classic Cox (1972) and Chevalier (1974) formalism used in 
each of the post-Larson (1974) papers listed above, highlighting an oversight which has 
lead, in some cases, to an overestimation of the gas ejection time ^gw- Following Cioffi, 
McKee k Bertschinger (1988), we discuss their alternative formalism for the evolution of 
the SNR thermal energy Stiiit), superseding the earlier Cox and Chevalier work with its 
more sophisticated treatment of radiative cooling in the hot SNR interior. Sections 3 and 
4 provide the relevant equations governing the evolution of eth(0 these competing 

formalisms. A direct comparison of the differences in some of the model predictions based 
upon the selection of thermal energy formalism is given in Section 5 and summarised in 
Section 6. Our aim is not to analyse each and every consequence of this selection, but 
only to highlight a few important ones, and in particular, to demonstrate the influence of 
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eth(0 upon ^GW predictions. This is sufficient for our purposes as many of tfie predictions 
(e.g. mass and metaiiicity of gas ejected) are related in a simple manner to ^gw- 



2. The Classical Wind Model of Ellipticals. For convenience we have adopted the clas- 
sical one-zone wind model of ellipticals as presented by Matteucci & Greggio (1986) and 
Matteucci & Tornambe (1987). A brief outline follows, with a more detailed account- 
ing of the relevant equations and input ingredients pertaining to our package (entitled 
MEGaW=Metallicity Evolution with Galactic Winds) to be presented elsewhere (Gib- 
son 1995), although the reader is directed to the above references, as well as those of 
Arimoto & Yoshii (1987) and Angeletti & Giannone (1990), for a complementary intro- 
duction to the subject. 

The gas mass Mg in the ISM is depleted by star formation and replenished through 
ejection from stars, and is governed by 

^ = -«0 + £(0, (1) 

where i^it) is the star formation rate and E{t) is the total gas ejection rate from all stars. 
Coupled to the gas mass equation 1, the mass of metals Mz = ZM^ in the ISM evolves as 

^^^ = -Z{mt) + Ez{t). (2) 

where Ez{i) is the total ejection rate of both processed and unprocessed metals. The 
relations for E{t) and Ez{i) are adapted from Greggio & Renzini (1983) and Matteucci 
k Greggio (1986). 

The nucleosynthesis prescriptions utilised are a function of both mass and metaiiicity. 
For the purpose of this paper we have used the yields of Arnett (1991) for Type II SNe, 
Renzini & Voli (1981) for single low and intermediate mass stars, and Thielemann, Nomoto 
& Yokoi (1986) for Type la SNe. Details related to the yields and the assumed SNe 
progenitors used in MEGaW can be found in Gibson (1995). 

The star formation rate i^it) is taken to depend upon the gas mass Mg(t) such that 

^(t) = i^M^(t) [Me/yr], (3) 

where the astration parameter z/ is the inverse of the star formation timescale, which for 
a galaxy of mass Mg is given by (Arimoto & Yoshii 1987) 

// = 8.6(Mg/1O^2^0)-°-^^' [Gyr-i]. (4) 

Recall that for t > ^gw, the star formation rate tjj(t) is taken to be zero. 

We use a singular power law of slope x = 1.35 (Salpeter 1955) or a; = 0.95 (Arimoto 
& Yoshii 1987) for the IMF, normalised to unity over the mass range mi = 0.1 Mq to 
m„ = 120.0 M0. Glisten & Mezger (1983) stellar lifetimes are used, and their respective 
remnant masses are taken from Maeder (1992). 
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As alluded to previously, the time at which global gas ejection occurs in the classical 
wind model ^gw is set by the residual thermal energy of the interstellar gas _Eth overcoming 
the binding energy of said gas i7g (Larson 1974), where 

^th(0= t eti,(t-t')Rst^(t')dt'. (5) 



t' is the SN explosion time. i?sN(^') is the total SNe rate (Type la and II) at time t' , and 
is taken taken from Greggio & Renzini (1983) and Matteucci & Tornambe (1987). Finally, 
Eth represents the evolution of the thermal energy content in the interior of a SNR. 

rig(i) is sensitive to the assumed dark matter distribution. Matteucci (1992) finds 
that for realistic diffuse halo distributions, the effect upon ^gw, and the chemical evolu- 
tion in general, is negligible. Hence we shall defer discussion of the diffuse dark matter 
halo model to her paper, as well as to recent work by Bressan, Chiosi & Fagotto (1994) 
and Forman, Jones & Tucker (1994). For our purposes, we shall use the binding energy 
formalism appropriate for a two-component spheroid (dark matter -|- luminous matter) 
in which the components are distributed similarly, and the ratio of the masses of dark to 
luminous matter is denoted R. Thus, using Saito (1979) and Matteucci (1992), we have 

^ . s ^ MM ) f MM ) \ 

where the initial binding energy of the system JIg of initial mass Mq = Miuin(l -|- R) (in 
grams) is 

1)g = 3.41 X lO-^M^-45 [erg]. (7) 

We can see from equations 6 and 7 that the ratio of dark to luminous mass R has the 
effect of increasing the binding energy of the gas by a factor of (1 -|- R)^'"^^ and the virial 
radius of the spheroid by (1 -|- i?)°'^^ (^-g- equation (6) of Matteucci 1992). What this 
means for predicted epochs of gas ejection will be addressed in Section 5. 

Finally, the exact formalism adopted for eth impacts greatly upon the predicted 
values of ^gw, and hence upon the predictions of mass of metals ejected at ^gw, as we 
shall discuss in Section 5. First though we shall introduce the competing formalisms 
available describing the evolution of eth, including the classic formulation of Cox (1972) 
and Chevalier (1974) (hereafter referred to as the "CC formalism"), as well as the newer 
work of Cioffi, McKee & Bertschinger (1988) (hereafter, the "CMB formalism"). It is 
important to stress that all of the previous wind models listed earlier (e.g. Saito 1979; 
Arimoto & Yoshii 1987; Matteucci & Tornambe 1987; Angeletti & Giannone 1990; Bressan, 
Chiosi & Fagotto 1994; Elbaz, Arnaud & Vangioni-Flam 1994; Gibson 1994) use the CC 
formalism. The work presented herein is the first to explore usage of the CMB form. 



4 



3. Supernova Remnant Thermal Energy Evolution: CMB Formalism. The evolution of 
a SNR can be characterised by three dynamical phases: (i) free expansion until the mass 
of interstellar material swept up reaches that of the SN ejecta, at which time a shock wave 
forms; (ii) adiabatic expansion which continues until the radiative cooling time of newly 
shocked gas equals the expansion time of the remnant; and, (iii) formation of a cold dense 
shell behind the front which begins when some sections of the shocked gas have radiated 
most of their thermal energy, being further compressed by the pressure of the remainder 
of the shocked material. Comprehensive reviews of the relevant shock dynamics include 
Ostriker k McKee (1988) and Draine k McKee (1993). 

Using a Id-Lagrangian hydrodynamics code, Cioffi, McKee k Bertschinger (1988) 
(CMB) followed the dynamical and thermal evolution of SNRs from the initial ejecta 
phase, into the adiabatic phase, and through to the phase in which radiative cooling of 
the hot interior becomes important. By insisting that the expansion of the shock radius 
Rs obey a power law in time (i.e. Rs oc t'^), CMB identified each of the phases simply by 
adopting the appropriate slope r] for each stage, the value being governed by a series of 
fits to their numerical models. 

The initial free expansion phase enters the adiabatic phase once the swept-up in- 
terstellar material begins to dominate the remnant mass, at which point the shock front 
slows down. The familiar Sedov- Taylor solution for a self-similar adiabatic shock with 
ri^r^ = 2/5 is encountered {t in s): 

i2.= (^)'V [cm], (8) 

where R^ is the radius of the outer edge of the SNR shock front, Eq = EqX 10^^ erg is the 
initial blast energy (i.e. Eq represents said energy in units of 10^^ erg), pQ is the ambient 
mass density (g/cm"^), and the constant ^ = 2.026 (appropriate for 7 = 5/3) (Ostriker k 
McKee 1988). 

Radiative cooling of the shocked material leads to the formation of a thin, dense 
shell at time ^f, where 

4f = 3.61 X [yr], (9) 

which "snowplows" through the ISM, driven from behind by the pressure of the hot SNR 
interior, as well as its own momentum. is a parameter which takes into account the 
infiuence of metallicity in determining the shell cooling time ((m = ^1^®)- 1^0 is the 
Hydrogen number density (cm~"^) in the surrounding ambient gas. A recent study of thin 
shell formation timescales corroborates the use of equation 9 for typical uniform ambient 
ISM densities (Franco et al. 1994). 

Before ^f, at time tp^s ~ 0.374f, radiative losses begin to affect the evolution of 
the remnant. At this point, the post-shock fiuid velocity approaches the shock velocity, 
entering the the so-called "pressure-driven snowplow" (PDS) phase. If one ignores cooling 
of the hot interior, the dynamics of the shell are governed by Rs oc ^2/7 (Ostriker k McKee 



5 



1988). CMB found a more suitable fit in the PDS stage using r]pds = 3/10, as opposed to 
this standard value of 2/7. Specifically, they recommend the PDS solution 



4t 

3^pds 



-1^^3/10 

3, 



where the radius at the beginning of the PDS stage is 

7? _ 14 nc2/'^„-3/'^/'-l/7 



[pc]. 
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Recall that tp^s in equation 10 is given by tp^s ~ 0.374f- 

Fitting to their simulations, CMB then present the following empirical representation 
for the evolution of the interior thermal energy of a SNR eth(0 times t < a24f , where 
the constant a2 = 1.169: 

0.717£;o(l - ai) 



1 



t \14/5' 



R, 
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1 



1/5 



1/9' 
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where ai = 0.398. Thus, the thermal energy evolution during the early dynamical phases 
of a SNR's history can be determined by coupling equations 8 (0 < t < ipDs), 9, and 10 
(^PDS < ^ < (i2tsf), with equation 12. 

The interior gas continues to lose energy by pushing the shell through the ISM 
(eth oc t~^l^), but there is an additional source of energy loss in the PDS phase, that 
being radiative cooling (eth t~'^^^). CMB provide the first analytical treatment of SNR 
evolution which incorporates explicitly these radiative cooling losses from the interior. The 
late-time behaviour of their remnants does approach the expected "momentum-conserving 
snowplow" (MCS) solution of Oort (1951), with Rg oc t^^'^ (at time t > ^mcs), eventually 
merging with the ISM when the internal driving pressure vanishes due to the radiative 
cooling in the interior, provided that this radiative cooling has not robbed the interior of 
the necessary hot gas to continue driving the shell through the MCS stage. In general 
though, CMB find that their remnants merge with the ISM before ever entering the MCS 
phase because the interior pressure drops below the ambient ISM pressure at times much 
less than ^mcs- 

Similar to equation 12, CMB present an empirical form for the SNR interior thermal 
energy evolution eth(0 during this later PDS phase, valid for times t > a24f: 



£th(t) 



0.717£;o(l - ai] 



10 



1/5 



1/9- 
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In the CMB formalism, expansion of the remnant is halted at a radius i?merge once 
the interior pressure in the SNR is reduced to that of the pressure of the ambient ISM, 
thereafter merging and becoming indistinguishable from the surrounding ISM, at time 
Emerge, where 



t 



merge 



57ipds£o 



5/49^10/49.15/49 



(14) 
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and 

-Kmerge — 3.7i?pcls£g Kg . (15) 

Beyond Emerge, the thermal evolution is no longer governed by the combination of 
the expansion term (eth t~^l^) and the radiative cooling term (eth t~'^^^) in equation 
13, but by only the latter. This late-time scenario, denoted Model B2 in our work, is 
computed by fixing Rg in equation 13 to be i?merge (equation 15), for all times t > Emerge- 
Typically though, the radiative cooling only proceeds for a time ^coob where 

Uoi ~ 550ipdsC/'', (16) 

until the hot interior gas has cooled to T ~ 10"* K, below which radiative losses become 
minimal (e.g. Sutherland & Dopita 1993; Saito 1979). Our Model B3 follows this pre- 
scription, depositing the residual SNR interior thermal energy at time ^cool into the ISM 
with no subsequent energy loss. 

While Model B3 is our preferred conservative model, alternative late time thermal 
energy evolution scenarios have been explored. Model Bi is similar to the model used 
by Larson (1974), in which radiative losses are negligible after the remnant has merged 
with the ISM (at time t = Emerge)- Model Bq is the opposite extreme; in this case, SNRs 
continue to expand in the PDS phase and cool radiatively for all times. As will be shown 
in Section 5, Model Bq leads to unfeasible results. 

Figure 1 highlights the differences in the late-time thermal evolution eth for the four 
different models explored using the CMB formalism. The example shown represents a 
SN blast of initial energy Eq = 10^^ erg, expanding into an ambient homogeneous ISM of 
number density no = 0.55 cm~"^, mass density po = 9.2x10"^^ g/cm"^, and solar metallicity. 
Model Bo is taken to expand and cool radiatively for all times (i.e. eth(0 ^ i~^'°'*). The 
remaining B models trace the evolution of Model Bq up until t = Emerge, the time at which 
pressure equilibrium between the interior and exterior ISM is reached, and expansion 
is halted. If radiative cooling can be neglected during the subsequent evolution, then 
eth(^) = eth(^merge)7 ^s was doue for Model Bi. In the more realistic scenario, radiative 
cooling reduces the SNR interior thermal energy such that eth(0 ^ until the time 

^cool at which point most of the hot gas within the remnant has cooled to T ~ 10"* K, at 
which point further decreases in the residual thermal energy are halted, as in Model B3. 
Model B2 continues to evolve like eth ^-0-44 fgj, times t > Emerge- 

Models Bi, B2, and B3 all stop expanding at time Emerge, once the internal pressure 
in the SNR is balanced by the ambient external pressure in the ISM. Following Larson 
(1974), we have also computed a parallel set of models (Models B^, B2, and Bg) which 
end their respective PDS phases at ^merge? where 

Cerge - 5.7 X W%'-'\r'i'-'-''Rh''- (17) 

^merge represents the characteristic time at which the SNRs collide and merge with other 
SNRs, thus halting further expansion. The star formation rate tp and galactic radius Rq 
are both set by the chemical evolution code MEGaW (Gibson 1995). The differences in 
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adopting equation 17 as opposed to equation 14 as the end of the PDS phase are usually 
minimal (see Section 5), but for completeness we do include this option. 

The dynamical (Rg (x t^ ) and thermal (eth (x f^) characteristics of the various models 
introduced above are summarised in Table I. 



4- Supernova Remnant Thermal Energy Evolution: CC Formalism. In the pioneering 
work of Cox (1972) and Chevalier (1974) (together referred to as CC, henceforth), as for 
the CMB formalism, the dynamics of the early stages of expansion are represented by the 
Sedov- Taylor similarity solution for an adiabatic blast wave (equation 8). 

In contrast with equation 9, the complete radiative cooling time of the SNR shell tc 
is given by 

= 5.7 X lO^e^/^'n-'/'' [yr]. (18) 

tc from equation 18 is typically a factor of two larger than 4f as given by equation 9 in 
the CMB formalism, a difference found to be negligible in predicting the epoch of wind 
ejection ^gw- 

The thermal energy in the interior of a SNR during the adiabatic expansion phase 
(t < tc) is simply given by 

eth(0 = 0.72^0, (19) 

where, as before, the initial blast energy Eq = Sq X 10^^ erg. 

Unlike that for the CMB formalism, CC ignore the transition from the Sedov- Taylor 
solution to the pressure-driven snowplow, which following the notation of CMB occurs at 
^pds < t < tsf (see Figure 1). Hence, tp^s and tc coincide in the CC formalism, and a 
discontinuity is encountered at the entrance to the PDS phase. 

The dynamics of the SNR shell during the subsequent PDS phase follows an Rg oc 
^0-31 law (Chevalier 1974). Neither Cox (1972) nor Chevalier (1974) have included radiative 
cooling of the hot interior, assuming that the cooling time is always longer than the 
dynamical time. In their formalism, the PDS (t > tc) SNR interior thermal energy evolves 
as 

eth(0 = 0.22^o(i/ic)"°-^'. (20) 

The thermal energy in the CMB formalism cools more rapidly (eth ^ as opposed to 

Sth i~°'^^) due to this inclusion of radiative cooling. Again, see Figure 1 for an example. 

As discussed in the previous subsection, we have chosen to test two different late- 
time thermal evolution histories. First, we allow all remnants to expand in the PDS phase 
for all times t > tc- This scenario is denoted Model Aq and is the procedure adopted 
by all previous researchers, including Arimoto & Yoshii (1987), Matteucci & Tornambe 
(1987), Angeletti & Giannone (1990), Bressan, Chiosi & Fagotto (1994), Elbaz, Arnaud 
& Vangioni-Flam (1994), and Gibson (1994). Again, a more appropriate approach is to 
halt the expansion at Emerge (Model Ai) or i^erge (Model A^ and Larson (1974)), as the 
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SNRs will not expand freely for all times; they will either come into pressure equilibrium 
with the external ISM, or start coming into contact with neighbouring shells. 

The power law behaviour of the dynamics of the shells and interior thermal energy 
for the models adopted under the CC formalism are listed in Table II. 



5. Discussion. In column 2 of Table III we list the epoch of global galactic winds ^gw for 
three different initial galactic masses Mg(0) based upon the assumption of the models for 
late-time evolution of the SNR interior thermal energy eth(0 described in Sections 3 and 
4. Recall that previous work (e.g. Arimoto & Yoshii 1987; Matteucci & Tornambe 1987; 
Angeletti & Giannone 1990; Bressan, Chiosi & Fagotto 1994; Elbaz, Arnaud & Vangioni- 
Flam 1994; Gibson 1994) adopted the CC formalism Model Aq. Columns 3, 4, and 5 show 
the gas fraction ejected at ^gw, the metallicity of said gas, and hence, the mass of metals 
ejected from the galaxy. The tabulated results are for spheroids without a dark matter 
component (i.e. i? = in equation 7). 

A comparison of the predictions of Models Aq and Ai proves particularly enlighten- 
ing, as they provide an extreme example of the sensitivity to eth(0- Both follow the CC 
formalism used in the above papers (i.e. Sthit > ^c) ^ but whereas Model Aq as- 

sumes that all SNRs expand and cool in the PDS phase ad infinitum, this expansion (and 
hence cooling) is constrained in Model Ai to halt at time Emerge once the shells come into 
pressure equlibrium with the ambient ISM. Alternatively, the expansion can be taken to 
halt once the SNRs start to overlap at time i^erge (Larson 1974). As can be seen in Table 
III, this scenario (Model A^) leads to ejection timescales of the same order of Model Ai, 
and thus for brevity, will be neglected henceforth. Regardless of initial mass, the epoch of 
gas ejection for Model Aq is five to six times later than that for Model Ai. For massive 
galaxies this implies that ~20 times the amount of gas and metals should be ejected at 
^GW in the Model Ai scenario. For the lower mass Ai models, the ejection time occurs 
before chemical evolution has enriched the gas beyond about solar metallicity. Thus the 
increase in metal ejection is simply due to the greatly increased mass of gas ejected. For 
the higher mass model, the chemical enrichment has reached approximately the same level 
as Model Aq. 

As might be expected from inspection of Figure 1, Model Aq is quantitatively similar 
to that of the CMB Model B2. As well. Model Bi parallels Model Ai in that expansion of 
the SNR shell halts at time Emerge (pressure equilibrium with the ambient ISM), although 
the resultant wind epochs occur somewhat later because the thermal energy content during 
the PDS phase is lower at a given time for the CMB form. 

One interesting result of using the CMB formalism is seen in the Bq models. These 
are similar to Aq models in that the SNR shells expand like Rg oc (i.e. eth R~^) ad 
infinitum. The key difference is that CMB also include radiative cooling of the hot SNR 
interior, the cooling term going as The infiuence this additional cooling has for the 

residual thermal energy available to heat the ISM gas at a given time is substantial. For 
masses > 10^° Mq, no galactic wind ever forms due to the increased cooling of the SNR 
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interior. This fact alone shows that a direct replacement of the older CC form for Model 
Ao with the newer CMB Model Bq is not valid. 

It is important to stress that all the models discussed here are at best a crude 
representation of the reality of ISM gas thermal evolution. In all the galactic wind papers 
discussed thus far, the ISM is assumed to be single smooth phase, whereas in reality 
the situation is much more complex. Three phase models (e.g. Ostriker & McKee 1988: 
cold cloud nuclei, surrounded by warm neutral gas, embedded in a hot thin corona) with 
substantial inhomogeneities are a step in the right direction, but difficult to quantify in a 
simple analytical form for use in galactic wind modeling (C.F. McKee 1994, priv comm). 
Hydrodynamical simulations of galactic winds from ellipticals, with (e.g. Missoulis 1994) 
and without (e.g. Ishimaru et al. 1994; Friaca & Terlevich 1994) cloudy multi-phase 
ISMs, are only now becoming feasible. A recent analytical multi-phase wind model has 
been proposed by Ferrini & Poggianti (1993), albeit at the expense of a large number of 
free parameters. 

Of even greater importance than the imposed single phase elliptical ISM may be the 
influence of magnetic field and cosmic ray energies in setting the timescale for gas ejection. 
It seems apparent that at least for the case of the local galactic disk (e.g. Duric 1988, and 
references therein), the magnitude of cosmic ray, thermal gas, and interstellar magnetic 
field energy densities are comparable to one another. It has been demonstrated that cosmic 
rays travelling through a magnetised plasma with mean velocity greater than the Alfven 
speed can become coupled to the plasma by the emission of magnetohydrodynamic waves, 
thereby "dragging" along the thermal gas as the rays escape the galaxy (Ipavich 1975). 
For galaxies of mass ~ 10^° Mq and magnetic field strengths similar to that found locally 
in our Galaxy (i.e. B ~ 5//G), Ipavich (1975), Zank (1989), and Fichtner et al. (1991), 
for example, have demonstrated that cosmic ray-driven galactic mass loss rates of 1-10 
M0/yr can be generated. Even if such winds are sustained for only 0.1-1 Gyr, the ejected 
mass of gas and metals would be substantial. 

Magnetic fields in the ISM also play another, secondary, role in SNe-driven wind 
models, by infiuencing the three dimensional dynamical evolution of an individual SN 
remnant. If the magnetic field strength is high enough for the outgoing shock to be weak, 
then most of the initial explosion energy Eq is used up doing work against the field and the 
initial expansion is halted rapidly, and no adiabatic Sedov- Taylor phase is encountered. 
Further expansion is then only possible along the direction of the interstellar field, with 
the late-time evolution being characterised by contraction of its equatorial regions and 
stretching in the axial direction (Insertis & Rees 1991). This "strong" field (e.g. B ~ IG) 
scenario is expected to be relevant for SNe situated in the centers of galaxies, and more 
importantly, for those in active galactic nuclei. When the interstellar Alfven speed is 
small compared to the velocity of the SN ejecta (i.e. the so-called "weak" field scenario), 
the initial phases are similar (i.e. the free expansion and Sedov- Taylor regimes) to the 
unmagnetised case. The late-time evolution is still governed by the halting of the equatorial 
expansion and unimpeded expansion of the shock along the field lines, again resulting in 
"barrel-shaped" remnants aligned with the magnetic field (e.g. Ferriere & Zweibel 1991; 
Tomisaka 1992). 
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Virtually all of the work just discussed has been tailored to studies of the Milky 
Way and other local spirals, simply because the dynamo- action necessary to generate 
large-scale galactic magnetic fields (Parker 1971) has only been confirmed in these sites. 
Recently though, depolarisation measurements of strong radio lobes within elliptical hosts 
(Garrington 1990) have implied magnetic field strengths of ~ 10//G, and the origin and am- 
plification of these fields within the framework of the turbulent mean-field dynamo theory 
demonstrated (Lesch & Bender 1990), although whether the properties of these confined 
fields can be extrapolated and assumed to permeate the ISM of "normal" ellipticals is still 
a matter of debate. 

The CC formalism outlined in Section 4 approximates the effects of weak (i.e. 
~ 3//G) magnetic fields by adding a term to the fiuid equations (Chevalier 1974). In 
Chevalier's analysis, the magnetic field was never energetically important unless the post- 
shock magnetic pressure exceeded the thermal pressure so that the magnetic field impeded 
shock compression. CMB neglect magnetic fields altogether but claim that they are im- 
portant for global evolution only if they inhibit the formation of the thin shell, which for 
the weak field scenario seems to be true qualitatively (Insertis & Rees 1991). 

The entire question of the role played by cosmic rays and magnetic fields in ellipticals 
is an unanswered one and beyond the scope of this paper; thus we refer the reader to the 
many studies described above for discussions of their importance for wind generation in 
spiral discs. We re-iterate then, that in determining the timescale for gas ejection ^gw, we 
have restricted ourselves to monitoring the infiuence of the gas thermal component only, 
with the caveat that substantial modification may be expected when a full accounting of 
cosmic ray and magnetic field energy is incorporated in future studies. For what follows, we 
provide the ejected masses at ^gw, and note the mass of gas and metals which accumulates 
in the ISM in the quiescent post-wind phase (i.e. ^gw < t < = 15.0 Gyr). We 
then speculate that some fraction of this reservoir may be ejected by one or more of the 
following: ram pressure stripping by other galaxies or the ICM; smaller late-time bursts 
of wind inherent to the galaxy; cosmic ray-driven steady mass-loss; post-^GW Type la 
SNe-driven steady mass-loss; etc. 

As alluded to in Section 3, of the options presented herein, we feel Models B3 or B3 
are currently the best compromise for modeling eth- Expansion of the shell is halted at the 
appropriate time (pressure equilibrium with the ambient ISM or contact with neighbouring 
shells), and radiative cooling proceeds until the temperature drops to ~ 10"* K, beyond 
which the cooling becomes inefficient. Such a form for the thermal energy evolution leads 
to wind epochs which are two to three times earlier (see Table III) than those found 
by previous researchers who have used Model Aq, and results in three to four times the 
mass of gas and metals being ejected from each elliptical. Putting this result into a full 
observational context is done in Gibson (1995), although two examples of the implications 
of adopting the CMB versus CC formalism are now discussed. 

One obvious application of galactic wind models is in predicting the contribution 
of wind ejecta to the mass of gas and metals (in particular. Iron) observed in clusters of 
galaxies (e.g. Matteucci & Vettolani 1988; Renzini et al. 1993). We will not attempt a 
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full analysis of the ICM abundance question here (see Gibson (1995) for more details), but 
will illustrate one example (Iron in the Virgo Cluster ICM) of the sensitivity of predicted 
ejected Iron masses to the assumed eth formalism. Our goal is to attempt to account for 
the observed Virgo ICM Iron mass (i.e. (1.6 - 2.0) X 10^° Mq - Arnaud 1994) within the 
confines of this SNe-driven galactic wind model. 

Following the procedure outlined in Elbaz, Arnaud & Vangioni-Flam (1994), we first 
generated a series of power law fits relating a spheroid's present-day V-band luminosity to 
the mass of ejected Iron. The former was estimated using the fuel consumption theorem 
(FCT) of Renzini & Buzzoni (1986), coupled with the typical observed rest frame (B-V)o 
colours from Brocato et al. (1990). As in Padovani & Matteucci (1993), a multiplicative 
correction factor of 1.4 was applied to the FCT-determined Ly. Two scenarios for the 
Iron ejecta were adopted. First, only the Iron ejected at ^gw was considered, as this 
"conservative" quantity is determined reasonably well by the wind model. Second, we 
estimate the maximum Iron available for deposition to the ICM (via one or more of the 
physical processes listed above) by adding in the mass of Iron which accumulates in the 
ISM reservoir during the post-^GW regime. With these Mp^ versus fits in hand, and 
the assumption of iv!E+so = 1-15 X 10^^ L® {Hq = 50 km/s/Mpc) (Elbaz, Arnaud k 
Vangioni-Flam 1994), we can predict the ICM Iron mass due to Virgo cluster ellipticals 
and lenticular s. 

Table IV shows the predicted Virgo ICM Iron masses for both the "conservative" 
(i.e. Mpg^) and "maximal" (i.e. Mpg^"*"**^) models. Parallel predictions are made for the 
classic SNR thermal energy evolution Model Aq and, for the first time, the newer CMB 
Models B3 and Bg. We have adopted both the Salpeter (1955) power law IMF of slope 
X = 1.35, and the fiatter x = 0.95 slope preferred by Arimoto & Yoshii (1987). 

Without resorting to the truncated IMF model of Elbaz, Arnaud & Vangioni-Flam 
(1994), the immediate conclusion to be drawn from Table IV is that the conservative 
X = 1.35 IMF models are unable to account for the (1.6 — 2.0) X 10^° Mq of Iron observed. 
Model B3 is certainly an improvement upon Model Aq, which underproduces Iron by a 
factor of ~ 27, but is still too low by a factor of ~ 6. The "maximal" Model B3 with IMF 
X = 1.35 is only a factor two lower than the observations. The fiatter Arimoto & Yoshii 
(1987) IMFs coupled to the CMB formalism Models B3 and B3 can easily reproduce the 
Virgo Iron abundance without resorting to any additional post-^GW mass loss. It should 
be noted that Model Aq with the fiatter IMF can only recover the observations provided 
substantial late-time mass loss is invoked. The recent work of Elbaz, Arnaud & Vangioni- 
Flam (1994) does not support post-^GW mass-loss at such a level. 

The above discussion is by no means exhaustive, and is meant simply to be illustra- 
tive. Again, a more complete analysis is forthcoming (Gibson 1995). 

Thus far we have ignored the effects of including a dark matter component dis- 
tributed similarly to the luminous component. Specifically, we have assumed the ratio of 
dark matter mass to luminous mass R in equation 7 to be zero. Taking R = b 10, 
Matteucci (1992) rightly demonstrated that use of the Model Aq formalism resulted in 
wind epochs of ~ 10 Gyr for initial galactic luminous masses of > 10^^ M© (e.g. her 
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Models CI and C2). Recalling that active star formation does not cease until after the 
global wind (Larson 1974), Matteucci (1992) noted that based upon the "late" wind times, 
active star formation (i.e. tp = 1 ^ 10 M0/yr) should be observed in low redshift (i.e. 
z < 0.5) ellipticals. Since this magnitude of star formation is not observed (e.g. Sandage 
1986; Burkert & Hensler 1989) Matteucci (1992) concluded that the dark matter cannot 
be distributed like the luminous matter, but in a diffuse halo. While our aim is not to 
argue strongly against this conclusion, we would like to illustrate at least what the impact 
of changing from Model Aq to Models Ai or B3 would have. 

Following Bressan, Chiosi & Fagotto (1994), let us take i? = 5 in equation 7. As 
expected from Matteucci's (1992) work, using the conventional Model Aq for a 10^^ Mq 
model elliptical leads to ^gw ~ 13.5 Gyr, and an implied star formation rate tp 4 
M0/yr at 2; ~ 0. Star formation at this level is what is observed currently in the local 
solar neighbourhood (i.e. late-type spiral disks), but not in typical nearby ellipticals. 
Using Model Ai yields ^gw ~ 0.8 Gyr, while Models B3 and B3 give ^gw ^ 2 Gyr. As 
we can see, using a SNR interior thermal energy evolution form which properly halts the 
expansion of the shell once pressure equilibrium between the interior and exterior ISM 
has been achieved (or once coming into contact with neighbouring shells), leads to wind 
epochs which occur early in the spheroid's lifetime, and hence active star formation at 
2; < 1 is thereby avoided. Thus, based upon this chemical evolution argument alone, 
one should not dismiss outright the "dark matter distributed similarly to the luminous 
matter" scenario, as done by Matteucci (1992). Indeed, the recent work of Berlin et al. 
(1994) shows that there does exist many ellipticals whose photometric and kinematical 
data are consistent with the model in which the dark and luminous matter are distributed 
similarly, or even with an absence of dark matter altogether - five of the six ellipticals in 
their small sample fall into this category. This does not invalidate the diffuse dark halo 
model, as stressed by both Saglia, Berlin & Stiavelli (1992) and Berlin et al. (1994). 
We have illustrated simply that the dynamically dominant dark matter scenario can be 
accommodated within the framework of the classic SNe-driven wind model, provided one 
account properly for the late-time evolution of the SNRs internal thermal energies. 



6. Summary. We have demonstrated the dangers in blindly adopting the SNR interior 
thermal energy evolutionary scenario in which each shell forever expands and cools in the 
pressure-driven snowplow phase. The most extreme example of this is seen when radiative 
cooling of the hot dilute interior is incorporated via the Cioffi, McKee & Bertschinger 
(1988) formalism; no galactic winds were found to develop in most of the model ellipticals 
unless the shell expansion was halted once coming into pressure equlibrium with the ambi- 
ent ISM. Conversely, very early galactic wind epochs were predicted when the Cox (1972) 
and Chevalier (1974) form was coupled to this pressure equilibrium constraint due to the 
absence of interior cooling. We suggest the compromise scenario (our Models B3 and B3) 
in which the late-time evolution of the thermal energy goes as eth(0 ^ (cooling and 

radiative cooling working together) until the shell either comes into pressure equilibrium 
with the ISM (Model B3) or contact with neighbouring shells (Model Bg). Cooling by 
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radiative processes alone then leads to eth(0 ^ ^-0-44 ^^til the interior cools to T ~ 10"* 
K, beyond which the efficiency of radiative cooling becomes negligibly small. Use of Model 
B3, as opposed to Model Aq which has been used predominantly in the past, leads to the 
onset of the galactic wind in approximately half the time. Consequently, the mass of gas 
and metals ejected is a factor of three to four times greater in the "newer" formalism. 

When a dark matter component distributed similar to the luminous component is 
introduced, use of the Model Aq formalism leads to late wind epochs (e.g. ^gw ~ 5 ^ 10 
Gyr) for typical mass-to-light ratios, implying active star formation in 2; < 1 ellipticals. 
Since star formation of this magnitude is not observed, indirect arguments have been made 
in the past that the dark matter in ellipticals must then be distributed in a diffuse halo. 
Applying Models Ai or B3 we find that the predicted wind times occur significantly earlier 
in the spheroid's evolution (e.g. ^gw ^ 4 Gyr, or 2; > 1), and thus this indirect argument 
alone is not sufficient evidence for a diffuse dark matter distribution in ellipticals. 

A conservative model for the Iron abundance in the Virgo Cluster ICM is proposed 
whereby a fiatter than Salpeter (x = 1.35) IMF of slope a; ~ 1, coupled with the CMB B3 
SNR thermal energy evolution models, can easily explain the measured metal mass. The 
classic Model Aq can be salvaged only if a substantial amount of post-^GW mass loss is 
invoked. The quantity of mass loss necessary seems at odds with the recent work of Elbaz, 
Arnaud & Vangioni-Flam (1994). 

Finally, it is important to stress that our goal herein was to highlight simply a po- 
tential weakness in the SNR interior thermal energy evolution scheme adopted in previous 
galactic wind models. Using the newer Cioffi, McKee & Bertschinger (1988) formalism in 
place of that of Cox (1972) and Chevalier (1974) is a step in the right direction, but one 
must recognise that there are still inherent weaknesses in the one-phase elliptical models 
employed: First, we must re-iterate that the galactic ISM is not particularly well-described 
by the one-zone model. Two, three, and four-phase ISMs with varying degrees of inhomo- 
geneities (Ostriker & McKee 1988) are certainly closer to reality, but difficult to quantify 
in a simple analytic form (although see Ferrini & Poggianti (1993) for a recent attempt 
to do so). Second, one must also recognise that there could be other significant contrib- 
utors to an elliptical's ISM thermal energy, besides SNe. Magnetic field and cosmic ray 
energies play a substantial role in driving galactic winds in spirals, as their magnitudes 
are comparable to the SNe thermal energy, but a direct extrapolation to ellipticals is more 
problematic and will require future re-examination. Bressan, Chiosi & Fagotto (1994) have 
speculated recently that mass loss from stellar winds, and not SNe, is effectively the only 
important factor in driving a galactic wind, although this conclusion has been questioned 
by Gibson (1994). 

These inherent weaknesses of the single phase, classic galactic wind model do not 
imply that analytical modeling of this sort is without merit. Even in this "fiawed" state, 
its successful deployment in explaining, for example, the mass-metallicity relationship in 
ellipticals (e.g. Larson 1974) and abundance ratios in the ICM of galaxy clusters (e.g. 
Matteucci & Gibson (1994), and references therein), demonstrates its usefulness. 

The aim of this paper is only to draw attention to the fact that the simplest analytical 
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galactic wind models presented in the literature, thus far, have adopted a potentially 
inappropriate formalism for the SNR thermal energy evolution. Again, we do not presume 
to imply that use of the Cioffi, McKee & Bertschinger (1988) form provides a full and 
complete description of the ISM gas thermal evolution (due to the absence of cosmic ray 
and magnetic field energies, as well as the enforced ISM homogeneity). Our hope is that 
by illustrating the sensitivity of model predictions to the adopted SNR thermal energy 
formalism, we can stimulate further work towards quantifying analytically the infiuence 
of the aforementioned neglected ISM energy sources. 

We have now begun the process of coupling our chemical evolution galactic wind 
code MEGaW to a full spectrophotometric package, the results of which will be presented 
elsewhere (Gibson 1995). 
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TABLES 



TABLE I 
cmb formalism: 
Power Law Time Dependent Coefficients 
Governing the Dynamical (t], where oc f^) 
AND Thermal (z/, where eth oc t") Evolution of the SNRs 
FOR THE Models Described in Section 2.2.1. 



Model 


t < tpDS 


^2^sf ^ t ^ ^merge 


^merge ^ ^ ^ ^cool 


^ ^ ^cool 








z/ 






ri 




Bo 


0.40 0.00 


0.30 


- 1.04 


0.30 


- 1.04 


0.30 - 


1.04 


Bi 


0.40 0.00 


0.30 


- 1.04 


0.00 


- 0.00 


0.00 - 


0.00 


B2 


0.40 0.00 


0.30 


- 1.04 


0.00 


- 0.44 


0.00 - 


0.44 


B3 


0.40 0.00 


0.30 


- 1.04 


0.00 


- 0.44 


0.00 - 


0.00 


Note: T] 


and V vary smoothly from t = 


^PDS to t 


= a24f- e. 


g. Figure 


1. 
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TABLES 



TABLE II 
cc formalism: 
Power Law Time Dependent Coefficients 
Governing the Dynamical (t], where oc f^) 
AND Thermal (z/, where eth oc t") Evolution of the SNRs 
FOR THE Models Described in Section 2.2.2. 



Model 


t < 


<i t <i ^nierge 


^ ^ ^merge 








T] V 


Ao 


0.40 0.00 


0.31 -0.62 


0.31 -0.62 


Ai 


0.40 0.00 


0.31 -0.62 


0.00 -0.00 
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TABLES 



TABLE III 

Galactic Wind Epochs ^gw (Gyr) 
Gas Fraction /g, Metallicity z, and 
Mass of Metals Ejected Mz (m®) at ^gw- 
'A' Models Use CC Formalism 
'B' Models Use CMB Formalism 



Mg(0) = 1.0 X IQQ M, 



Model 




/g(^Gw) 


Z(iGw) 




Ao 


0.182 


0.068 


0.069 


4.74(6) 


A 


U.Uoo 


0.565 


0.012 


D.OOl^D ) 


A'l 


0.033 


0.565 


0.012 


6.83(6) 


Bo 


1.468 


0.003 


0.065 


1.72(5) 


Bi 


0.059 


0.374 


0.023 


8.68(6) 


B2 


0.229 


0.040 


0.080 


3.22(6) 


B3 


0.102 


0.197 


0.041 


8.09(6) 




Mg(o: 


) = 5.0 X 10^° M® 




Model 




/g(^Gw) 


Z(iGw) 


Mz(tGw) 


Ao 


0.570 


0.016 


0.077 


6.09(7) 


Ai 


0.118 


0.299 


0.030 


4.55(8) 


A'l 


0.104 


0.342 


0.027 


4.55(8) 


Bo 


n/a 


n/a 


n/a 


n/a 


Bi 


0.189 


0.157 


0.049 


3.81(8) 


B2 


0.615 


0.014 


0.074 


5.15(7) 


B3 


0.305 


0.062 


0.071 


2.20(8) 




Mg(o: 


) = 2.0 X 10^2 




Model 


^GW 


/g(^Gw) 


Z(iGw) 


Mz(tGw) 


Ao 


2.189 


0.005 


0.068 


6.47(8) 


Ai 


0.374 


0.103 


0.060 


1.23(10) 


A'l 


0.307 


0.147 


0.051 


1.50(10) 


Bo 


n/a 


n/a 


n/a 


n/a 


Bi 


0.477 


0.063 


0.070 


8.79(9) 


B2 


1.859 


0.006 


0.068 


7.88(8) 


B3 


0.813 


0.020 


0.077 


3.05(9) 
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TABLES 



TABLE IV 

Predicted Virgo Cluster ICM Iron Mass (in m®) 
Contributions from either Global Gas Ejection 

AT ^GW ALONE (Mp«w) OR ^gw-COMPONENT PLUS 

Complete post-^gw Ejection of ISM (Mp«w+*«). 



Eth Model 


IMF X 




/l/f*GW+*G 


Ao 


1.35 


5.93 X 10® 


4.28 X 10^ 


Ai 


0.95 


3.66 X 10^ 


2.64 X 10^° 


B3 


1.35 


1.47 X 10^ 


7.36 X 10^ 


B3 


0.95 


2.02 X 10^° 


7.10 X 10^° 




1.35 


2.52 X 10^ 


8.40 X 10^ 


B^, 


0.95 


2.50 X 10^° 


7.93 X 10^° 


Observed" 




(1.6 - 2.0 


) X 10^° 



Arnaud (1994) 



21 



FIGURE CAPTIONS 



Figure 1: Evolution of the interior thermal energy eth of a supernova remnant, normalised 
to 10^° erg. The sample models shown represent a total initial blast energy of 10^^ erg, 
expanding into an ambient interstellar medium of number density no = 0.55 cm~"^, mass 
density po = 9.2 X 10"^^ g/cm"^, and metallicity Z = 0.02. 
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